NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


NOISE  CANCELLATION' USING 
ADAPTIVE  ARRAYS 

by 

Constantinos  G.  fanikas 
June  1984 


Thesis  Advisor 


H.  Titus 


Approved  for  public  release;  distribution  unlimited 


T221536 


S  E  C  U  R I  ~  ^    CL13S!c'CA-|j>. 


=  A3E    ?Wh*n    Dxia   Entered; 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 


REP3R  ~  NO  M2  E  = 


2.  GCVT  ACCESSION  NO.  3   =>  EC  '  =>  i  E  N  _' 3  CATALOG  n  U  m  3  ?  R 


4.      T|  TLE  ("and  Subtitle; 

Noise  Cancellation  Using  Adaptive 
Arrays 


S.     TYPE   Or    SEP09T   4   =  ERiOC    COVERED 

Master's  Thesis; 
June  1984 


6       P£RPCPm:n5  0R3  JEPrP"S-M5E' 


7.  AuT-CR^J> 


Constantinos  G.  Manikas 


8       CONTRA"   0-=    3R»».T    N_w3ER's, 


9.     PERFORMING   :»G»v:i":/<    NAME    ANO    AOCPESS 

Naval  Postgraduate  School 
Monterev,  California  93943 


10.      PROG  R  Ai 


MENT  ==oj  £•: 


ARE*    1    "VOR<    JNlT    >   JMBER5 


Naval    Postgraduate    School 
Monterey,     California    93943 

12       REPO-fT  I'Ts 

June    1984 

i    1 3 .      NUMBER  OP   PAGES 

I                TO 

r.:-r.   Controlling  Ollic*)      i    '5.      SECo^  ITY    CLASS 

Unclassi  tied 


I5a      DECLaSSi^    Ca"'ON     CO*' 


!     5. 


Approved  for  public  release;  distribution  unlimited 


"!       0\%~  a    3'JTICN    "i'Ev".-.*  ~<r    eoe'rec:    errered  <r.   Block   20,    II  dUferent   from   Report; 


IS.      SUPPL  EMEN  T  A»V   NC~ES 


'9       KEY    *GR0S    '  Con'  rnue   '>r   revert')    tine    it  ne'eeeerv    rtntl   Identify   by    block   number) 

Adaptive  Noise  Cancelling 

Adaptive  Arrays 

Hand  Contained  Beamf ormers 


*3S"=;:' 


\erv    sn.l    Itionltty    by    blorX    number* 


All  known  adaptive  be am formers  utilize  some  form  of 
automatic  minimization  of  the  mean  square  error.   High 
adaptation  rates  though,  exhibit  a  signal  cancellation 
phenomenon  leading  to  self-jamming  by  the  adaptive  antennas. 
This  effect  results  from  adaptive  interaction  between  signal 
and  interference  (i.e.,  jammer)  inputs  simultaneously 


^0     |    JAN    73     1473  EDITION   OP    ■    MOV   55   IS  OBSOLETE 

S    N    0)02-  LF-  0!  J-  56C  '  ]_ 


SECURITY    CLASSIFICATION    O^    ~S|S   PASS   f»n»n    Dmtm    Entered) 


SECURITY   CLASSIFICATION   OF   THIS   PACE  fWhmn  Datm  Enfrvd) 


20.   (Continued) 

received  by  an  adaptive  antenna.   This  research  investigates 
various  existing  ways  of  adaptive  beamf orming  for  noise 
cancelling,  and  signal  enhancement  from  simple  Adaptive 
Noise  Cancellers  to  Hard  Constraint  Adaotive  Beamf ormers . 


N   0102-  LF-  0)4-  6601 


SECURITY   CLASSIFICATION   OF   THIS  P  AGEfWhan  Data  Entarmd) 


Approved  for  public  release;  distribution  unlimited 


Noise  Cancellation  Using 
Adaptive  Arrays 


by 


Constantinos  G.  Manikas 
Lieutenant,  Hellenic  Navy 
Hellenic  Naval  Academy,  1974 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
June  1934 


ABSTRACT 

All  known  adaptive  beam  formers  utilize  some  form  of 
automatic  minimization  of  the  mean  square  error.   High 
adaptation  rates  though,  exhibit  a  signal  cancellation 
phenomenon  leading  to  self-jamming  oy  the  adaptive  antenna.3. 
This  effect  results  from  adaptive  interaction  between  signal 
and  interference  (i.e.,  jammer)  inputs  simultaneously 
received  by  an  adaptive  antenna.   This  research  investigates 
various  existing  ways  of  adaptive  beamforming  for  noisc-i 
cancelling,  and  signal  enhancement  from  simple  Adaptive 
Noise  Cancellers  to  Hard  Constraint  Adaotive  Beamformers. 
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I.    INTRODUCTION 


Adaptive  antennas  have  been  under  development  in  various 
forms  during  the  past  two  decades,  having  thus  far  proven 
themselves  capaole  of  rejecting  various  jamming  signals. 
Most  high  performance  radar  and  communication  systems  being 
designed  to  work  in  jamming  environments  currently 
incorporate  various  forms  of  adaptive  antennas.   The 
concurrent  development  of  spread  spectrum  techniques  >;i.  th 
the  adaptive  antennas  provides  a  set  of  technologies  'for  jam 
resistant  systems.   Both  of  the  above  technologies  are 
compatible  and  frequently  are  used  in  the  same  system.   The 
adaptive  antenna  is  to  attenuate  the  strong  jamming  signals 
as  they  appear  at  the  front  end  of  the  receiver;  then  spread 
spectrum  techniques  are  used  to  neutralize  a  large  number  Df 
weak  signals  that  may  not  be  eliminated  by  the  adaptive 
antenna,  and  to  recover  the  desired  signals. 

The  way  an  adaptive  antenna  works  is  to  pass  the  signal 
contaminated  with  noise  through  a  filter  that  tends  to 
suppress  the  noise  while  leaving  the  signal  unchanged. 
These  kind  of  filters,  known  as  adaptive  filters,  have  the 
ability  to  adjust  their  own  parameters  automatically  and 
their  design  requires  little  or  no  prior  knowledge  of  the 
signal  or  noise  characteristics. 


Adaptive  noise  cancelling  is  a  variation  of  optimal 
filtering  that  uses  an  auxiliary,  or  reference  input.   This 
input  is  then  subtracted  from  the  primary  input,  which  is 
composed  of  both  the  desired  signal  and  the  undesired  noise. 
This  procedure  of  filtering  and  subtraction  is  controlled  by 
an  adaptive  process,  giving  noise  reduction  which,  depending 
on  the  particular  application,  can  result  in  significant: 
noise  cancellation. 

The  design  of  an  adaptive  antenna  has  to  retain  signal 
preservation.   Unfortunately  some  adaptive  beam  formers  i  ". 
use  do  not  perform  well  in  certain  environments  where 
cor  related- s ignal  conditions  lead  to  partial  or  total 
cancellation  of  the  desired  signal  within  the  antenna. 

In  that  way,  although  our  primary  intention  is  the 
cancellation  of  all  the  undesired  interferences,  we  create  a 
self-jamming  of  the  adaptive  arrays  by  the  so-called 
phenomenon  of  signal  cancellation  resulting  from  the 
adaptive  interaction  between  the  signal  and  the  jammer  which 
simultaneously  are  received  by  the  adaptive  antenna's 
elements . 
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II .        ADAPTIVE    NOISE    CANCELLER 

A.        3 ASIC    CONCEPT 

The    oasic    block    configuration    of    an    Adaptive    No i 
Canceller    is    shown    in    Figure    2.1.       A    signal,    s,     is 


noise,    n    ,    uncor related    with    the    signal.       V'cn 


to  mea 


sign  r.    -  n  o  i  s e    ror.ii    t n e    p  r  i  :n  a  r  y    i  n  p  u 


i_     <_  <~> 


.  n  a 


i—  ."5  n  r"  — 


second    sensor    receives    a    noise,    n     ,     also    uncor related    w 

the    signal,    out    correlated     in    some    unknown    way    with    the 

noise,    n    .       This    second    sensor    provides    the    so-called 
'       o 

reference    input    of    tne    canceller.       Noise,    n ,  ,     is    filter: 
produce    tne    output,    y,    that    is    as    close    as    possioie    2 
replica    of    n    .       Output,    y,     is    now    subtracted    from    tne 
primary    input    to    give    the    final    system's    output. 


z    =    s    +    n      _    y 


(2.1) 


The  reference  input,  n,,  is  processed  by  an  adaptive 
filter  which  automatically  adjusts  its  own  impulse  responsi 
This  adjustment  is  accomplished  through  an  algorithm  tnat 
responds  to  an  error  signal  dependent  on  the  filter's 
output.   In  this  way,  using  tne  proper  adaptive  algorithm, 
the  filter  can  operate  under  changing  conditions  which 
minimize  the  error  signal. 
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The    error    signal,     e,    is    feeding    the    system's    output    sack 

to    t-he    adaptive    filter    in    order    to    ad]  ust    the    filter    through 
a    LMS     (Least    Mean    Squares)     adaptive    algorithm    and    minimizes 
the    total    system    output    power.       In    other    words,     in    an 
adaptive    noise    canceller    the    output    forms    an    error    signal 
for    the    adaptive    process. 

3.        ANALYTICAL    MODEL  — NOTCH    FILTER    RESPONSE 

In    order    to    analyze    the    adaptive    filter,     je    cho     = 
pure    sinusoidal     for    reference    input 


x .      =    C    cos(j    t    +     4> ) 

K  O 


(2.2 


Glover  [Ref.  1]  has  shown  that  under  this  condition  v"_ 
dashed  box  of  Figure  2.2(a)  can  be  approximated  by  Figure 
2.2(b)  . 

The  conventional  way  of  eliminating  such  sinusoidal 
interferences  is  through  the  use  of  a  notch  filter.   Figure 
2.3  snows  a  single  frequency  noise  canceller,  with  two 
adaptive  weights,  which  are  updated  through  the  equations: 


W1J+1  =  W1J  +  2UEJ  X1J 


(2.  3 


W2J+1  =  W2J  +  2ueJ  X2J 


(  2  .  4  ) 


The  sampled  reference  inputs  are 
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xLJ    =    C    cos ( w0]T    +    $) 


(2.5) 


x2J    =    C    sin(u>0JT    +    D) 


(2.6 ; 


Widrow  [Ref.  2]  and  [Ref.  3]  has  shown  that  the  single 
frequency  noise  canceller  has  the  properties  of  a  notch 
filter1  at  the  reference  frequency,  a  ,  having  the  ;:ans:- 
f unction  : 


Y  (2)  /X  (z)  =  H  (z) 


w  n  e  r  e 


H  (z! 


z2  -  2z  cos(2nj  n  L)  +  1 

o 


2  2 

z   -  2  ( 1  -  ,iC  )  2  cos  ( 2 


-1 


+  1  -  2u< 


The  zeros  of  the  above  transfer  function  are  lo< 
the  z-plane  at 


±J2tTO)   Q 

Z  =  e      o 


-1 


(2.8) 


and  are  inside  the  unit  circle. 

The  poles,  also  inside  the  unit  circle,  are  located  at 


z  =  (1  -  ye2)  cos  (2ir.jj  0  l)     ± 

o 


j[(l    -    2uc2)     -     (1    -    uc2)  cos2  (2ttojo^    l)]l/2 


(2.9) 


1 


Detailed  explanation  is  given  in  Appendix  A 
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Figure  2.4  indicates  the  location  of  poles  and  zeros  in 

the  z-plane.   Since  the  zeros  lie  inside  the  unit  circle, 

the  depth  of  the  notch  in  the  transfer  function  is  infinite 

at  the  interference  frequency,  u  .   The  sharpness  of  the 

o 

notch  is  determined  by  the  closeness  of  the  poles  to  the 
zeros  . 

Figures  2.5,  2.6  and  2.7  are  magnitude  plots  of  r_.ie 
transfer  function  Hfz)  using  as  C  the  value  of  one, 
different  ralues  of  tne  adaptation  constant  -  . 

C.   MOISE  CANCELLATION  IN  THE  NARROWBAND  CASE 

Consider  tne  case  in  which  ootn  signal  ana  jammer  [i.e. 

interference)  are  pure  sinusoids  of  known  frequencies,  :  ~, 

f_,  wnere 
J 


f  4    f 

"S  '    J  " 


Figure  2.3  shows  the  frequency  spectra  of  both  noise  and 
jammer,  and  Figure  2.9  indicates  the  same  spectra  after 
being  filtered  by  a  notch  filter  with  cutoff  frequency  the 
same  as  the  one  of  the  jammer,  f 

Comparing  the  two  preceding  Figures  2.3  and  2.9,  we  see 
how  effective  the  notch  filter  is  in  eliminating  sinusoidal 
interferences.   Of  course,  pure  sinusoidal  interferences  are 
seldom  found  in  nature,  and  still  the  main  disadvantage  is 
that  noise  and  jammer  have  to  be  of  different  frequencies. 
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Figure  2.4.   Pole-Zero  Configuration 
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III.   FROST  BASED  ADAPTIVE  3EAMF0RMER 

A.   CONSTRAINED  LEAST  MEAN  SQUARES  ALGORITHM 

Frost  [Ref.  4]  proposed  an  algorithm,  the  L.M.S. 
algorithm,  for  adjusting  an  array  of  sensors  in  real  time  to 
respond  to  a  desired  signal  while  3 iscr iminat ing  against 
noises.   Figure  3.1  illustrates  the  basic  set-up  for 
implementing  an  L.M.S.  algorithm.   It  is  assumed  t.nat 
desired  signal  arrives  in  plane  waves  from  a  chosen 
direction  called  the  look  direction. 

The  algorithm  iteratively  adapts  weights  in  :>rder  :o 
minimize  noise  power  at  the  array  output  while  maintai  .:  -j  .. 
chosen  frequency  response  in  the  look  direction.   Require- 
ments for  the  algorithm  is  a  priori  knowledge  of  the 
direction  of  arrival  and  the  frequency  band  of  interest. 
During  the  adaptive  process,  the  algorithm  progressively 
learns  the  statistics  of  noise  arriving  from  all  directions 
except  the  look  direction. 

A  major  advantage  of  the  constrained  L.M.S.  algorithm  is 
that  it  has  a  self-correcting  feature,  permitting  it  to 
operate  for  arbitrarily  long  periods  of  time  in  a  digital 
computer  implementation  without  deviating  from  its 
constraints  oecause  of  cumulative  roundoff  or  truncation 
errors  . 
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Noise  arriving  from  the  look  direction  may  be  filtered 
out  by  a  -suitable  choice  of  the  frequency  response 
characteristics  in  that  direction  or  by  some  external 
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B.   BASIC  PRINCIPLE  OF  THE  CONSTRAINTS 

L.M.S.  Algorithm  maintains  a  chosen  frequency  response 
in  the  look  direction  while  minimizing  output  noise  po .;  =  ■: 
because  of  a  simple  relation  between  the  loo:<  direction 
frequency  response  and  the  weights  in  the  array  of  Fig  ;r-- 
3.1.   Assuming  that  the  look  direction  is  chosen 
perpendicular  to  the  line  of  sensors,  identical  signal  . 
components,  arriving  on  a  plane  wavefront  parallel  10  the 

line  of  sensors,  appear  at  the  first  taps  simultaneously  and 

2 

parade  in  parallel  down  tne  tapped  delay  lines   nollcwirvj 

each  sensor.   However,  noise  waveforms  arriving  iron  other 
than  the  look  direction  will  not,  generally,  produce  equal 
voltage  components  on  any  given  vertical  column  of  taps. 
The  voltages,  signal  plus  noise,  at  each  tap  are  multiplied 
by  the  tap  weights  and  added  to  form  the  array  output. 
Thus,  as  far  as  the  signal  is  concerned,  the  array  processer 
is  equivalent  to  a  single  tapped  delay  line  in  which  each 
weight  is  equal  to  the  sum  of  the  weights  in  the 
corresponding  vertical  column  of  the  processor,  as  indicated 
in  Figure  3.2.   These  summation  weights  in  the  equivalent 


'Appendix  B  describes  the  tapped  delay  line  filter. 
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tapped  delay  line  must  be  selected  so  as  to  give  the  desirei 
frequency  response  characteristics  in  the  look  direction. 

If  the  look  direction  is  chosen  to  be  other  than  the  one 
perpendicular  to  the  line  of  sensors,  -hen  the  array  can  oe 
steered  either  mechanically,  or  electrically,  by  the 
addition  of  steering  time  delays  placed  immediately 
eacn  sensor. 

Mow  consider  the  array  problem  of  Figure  3.2. 

In  vector  notation  the  sample  tap  voltages  3.  re 
iden 1 1  £  led  as 


=  --_•- 


XT(])  =  [x]_  (j)  ,  x2  (j)  ,  .  .  .  ,  xkJ(j) 


and  the  sampled  weights  as 


T  .   A 
W  (J)  =  [wx(]),  w2(j),  ...,  wkJ(j) 


(3.2) 


where  j  represents  the  sample  number.   The  filter  output 
are  summed  to  form  the  beamformer  output 


y(j)  =  wT(j)  •  x(j)  =  xT(])  •  w(j 


(3.3) 


If  we  denote  by  C  the  constraint  matrix,  defined  as 


V-  LC-1,     C^,      ...,     C-r 


(3.4) 
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a  KJ  :<J  dimensional  matrix,  that  conveys  the  beam  f  orrn-er 
structural  info-rmat ion  needed  to  implement  the  constraints, 
then  the  constrained  L.M.S.  problem  may  be  stated  as 


Minimize  E[y2(j)]  =  E [ WT • X (j ) XT ( j ) • W] 


(3.5 


suoject    to      C    A    -  o<_ 


(3 


where    >•<_.   is    the    look-direction    response    -/ector    defined 


■£-     =     [fl'     r2 rJ 


(3.7) 


The  weights  are  updated  using  the  constrained  L.M.S. 
algorithm,  which  for  the  case  of  Figure  3.2  (assuming  two 
sensors) ,  can  be  written  as 

w1(j+l)  =  (1/2) [wx(j)  -  uy(j) xx(j) ] 

-  (l/2)[w2(j)  -  fiy('j)  x2(j)  ]  +  (fx/2)    (3.3) 

w2(j+l)  =  (l/2)[w2(j)  -  uy(j)x2(j)] 

-  (l/2)[w1(j)  -  uy(j)x1(j)]  +  (fx/2)     (3.9) 
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w3(j+l)  =  (l/2)[w3(j)  -  uy(j)x3(j)] 

-  (l/2)[w4(j)  -  uy(j")x4(j)]  +  (f2/2)   (3.10: 

w4  (j+1)  =  (1/2)  [w4(j)  -  uy(j)  x4 (]) 

-  (l/2)[w3(j)  -  uy(j)x3(j)]  +  (f/2)    (3.11 

The  per  f  o  cmance  that  a  FROST  oased  adaptive  beamformer 
exhibits  is  unity  gain  and  flat  frequency  response  in  the 
look  direction.   A  ieep  null  must  be  placed  at  the  jammer's 
frequency.   Figure  3.3  is  a  plot  of  the  signal  direction 
frequency  response,  and  Figure  3.4  is  a  plot  of  the  jammer- 
direction  frequency  response.   Figure  3.5  indicates  the 
beampattern  of  the  FROST  beamformer,  confirming  that  a  deep 
null  has  been  placed  on  the  jammer. 

C.   NOISE  CANCELLATION  IN  THE  WIDEBAND  CASE 

The  first  case  for  -examination  is  the  one  of  a 
sinusoidal  jammer  of  frequency,  f  ,  and  a  broadband  signal 


J 


of  center  frequency,  f  ,  where 


fa   =   fs  (3.12) 
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Figure  3.5.   FROST  A.B.F.  Beampattern 
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Figure  3.6  represents  the  spectrum  of  the  signal,  since 
Figure  3.7  indicates  the  spectrum  of  the  stronger  sinusoidal 
jammer.   Figures  3.3,  3.9  and  3.10  represent:  results  of 
passing  the  above  described  signal  and  noise  through  notch 
filters,  each  one  employing  a  different  adaptation  constant. 

Comparing  the  preceding  three  figures  we  see  how 
important  the  adaptation  constant,  u,  is  and  now  it  car.  ieai 
to  signal  aliasing  by  choosing  it  to  have  bigger  than  proper 
mag  ni tude . 
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IV.  DUVALL  BASED  ADAPTIVE  BEAMFQRMER 

A.   BASIC  CONFIGURATION 

We  have  already  illustrated  some  relatively  simple 
signal/jammer  scenarios  and  have  seen  thai:  not  only  the 
Adaptive  Noise  Canceller  (A.N.C.)  ,  but  even  the  most  riJ.4. 
constrained  beam  formers  (Frost  A.B.F.)  can  fail  to  preserve 
the  desired  signal  from  distortion. 

In  this  chapter  we  will  review  and  describe  another 
solution  to  the  noise  cancellation  problem  and  see  r.z .;    we 
can  preserve  the  desired  signal  at  the  price  of  some 
increase  in  beam  former  complexity. 

The  basic  configuration  of  the  DUVALL  based  beamf  ,■::-". 
[Ref.  5]  is  illustrated  in  Figure  4.1. 

Two  observations  were  made  by  Duvall  useful  in 
developing  the  Composite  Beamformer  (C.B.F.)  .   The  first 
observation  was  that  interaction  between  the  desired  signal 
and  the  ]amraer  is  the  root  for  the  signal  cancellation 
phenomenon.   Frost  experimentally  and  analytically  has  shown 
that  the  presence  of  both  signal  and  jammer  energy  is  a 
prerequisite  for  signal  cancellation,  and  the  output  signal 
containing  correlated  signal  and  noise  waveforms  is  the 
phenomenon  that  makes  the  signal  distortion  occur. 

The  second  observation  was  that  the  signal  plays  no  role 
in  the  optimum  (i.e.,  Wiener)  solution  calculation  in  a 
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perfectly  steered  Frost  A.3.F..   In  other  words,  the  look- 
direction  response  is  determined  exclusively  by  the  hard- 
constraints  and  not  by  the  characteristics  of  the  desired 
signal  . 

These  two  above  observations  were  addressed  for 
exclusion  of  the  signal  from  the  beamformer. 

Reconsidering  the  whole  structure  of  the  C.B.F.,  we  se- 
that  the  key  elements  are: 

1.  An  Augmented  Array.   The  array  has  been  augmente-i  sc 
that  su'oarrays  consisting  of  multiple  elements  appea 
in  place  of  the  individual  elements  of  the  :r  l j . 
array . 

2.  A  Preprocessor.   A  preprocessor  operates  upon  che 
received  signals  from  the  augmented  array  to  je-.  ^r ~  _■ 
an  environment  that  is  free  of  desired-signal  conten 


3.  An  Adaptive  Beamformer.   The  adaptive  beamformer  is 
required  that  can  be  constrained  to  control  che  look- 
direction  response,  while  nulling  all  jamming  signals. 

4.  A  Slaved  Non-Adaptive  Beamformer.   The  adaptive 
beamformer  derives  weights  that  are  copied  to  the 
slave  beamformer,  which  has  the  same  structure  as  the 
adaptive  beamformer  but  is  connected  directly  to 
selected  antenna  elements  and  is  used  to  implement  the 
computed  solution  and  recover  the  desired  signal. 

Three  of  the  above  key  elements,  the  array,  the 
preprocessor,  and  the  adaptive  beamformer,  afford  consider- 
able flexibility.   A  variety  of  specific  realizations  are 
possible.   The  slaved  beamformer  design  is  inflexible  in  the 
sense  that  it  mirrors  the  design  of  the  adaptive  beamformer. 


B.   PRINCIPLES  OF  THE  DUVALL  C.B.F. 

The  general  C.B.F.  structure  of  Figure  4.1  can  be 
specialized  to  yield  the  C.B.F.  shown  in  Figure  4.2.   Here 
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the  desired  signal  is  assumed  to  be  incident  from  oroadside 
so  we  can  neglect  beam  steering.   The  preprocessor  is 
realized  by  using  two-element  suoarrays  in  a  simple  element 
differencing  scheme.   The  snaring  of  elements  between 
subarrays  provides  very  efficient  use  of  elements;  only  a 
simple  element  is  needed  oeyond  those  ordinarily  required 
for  a  comparable  array. 

In  order  to  prevent  the  desired  signal  from  distortion 
(i.e.,  signal  cancellation),  the  signal  which  arrives  i  r  :  ■  r. 
tne  Look-direction  is  excluded  from  the  beamformer  /nere  the 
adaptive  process  takes  place.   Under  the  proposed  scheme, 
the  adaptive  process  is  only  used  to  determine  a  set  of 
weights.   These  weights  are  copied  into  a  separate, 
identical  processor  which  is  used  to  form  the  output  signal . 
The  C.B.F.  scheme  can  be  in  cooperation  with  any  existing 
adaptive  algorithm  including  the  Frost  one,  as  Figure  4.2 
addresses  . 

Due  to  the  subtractive  preprocessing,  the  look-direction 
signal  does  not  appear  at  the  Frost  adaptive  beam  former 
inputs,  whose  inputs  are  only  composed  of  the  jamming 
signal.   The  task  of  the  Frost  Algorithm  is  to  place  a  null 
at  the  gammer  direction  through  the  determination  of  an 
appropriate  set  of  weights  and  the  signal  has  now  no  effect 
at  all  on  the  determination  of  those  weights.   By  copying 
the  weights  the  slaved  processor  places  the  main  beam  of  the 
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antenna    pattern    to    the    constraints    established   by    the    Look- 
direction    and    verifies    that    a    null    has    already   been    placed 
in    the    jammer's    direction. 

Figure    4.3    represents    the    Duvall    C.B.F.     in    terms    of 
pnasor    notation.       The    jammer    components    received    by    the 
antenna    are    indicated    by    a    see    of    equal    amplitude    jniformly- 
spaced    phasors,    JQ/    J,,    J9,    J -,    and    J4-       The    phasor     inputs    :c 
the    Frost    b earn  former    are    J,     -    JQ ,    J„    -    J     ,     j       -    J         and 
j      -    j-,.       :c\ay    are    also    uniform-amplitude,    equally-spaced, 
and    separated    oy    tne    same    angles    as    the    received    jammer 
components,     JQ,    J,,     J2,    JU    and    J.. 

Since    the    relative    phase    angles    are    the    same    in    the 
slaved    processor    as    in    the    Frost    processor,     tne   beam 
pattern    notch    is    formed    at    the    proper    bearing    angle. 

The    pnasor    argument    applies    to    a    single    jammer    at    a 
single    frequency.       Linearity    and    superposition    apply    and 
show    that    phase    relations    are    preserved     for    multiple    jammers 
and    for    broadband    as    well    as    narrowband    signals. 

The    uniform    linear    array    provides    an    attractive 
structure    for    the    C.B.F.    because    there    is    the    option    of 
element    sharing    between    subarrays.       A    regular    array 
structure    is    not,    however,    a    prerequisite    for    the    C.B.F.. 
The    fundamental    requirement    is    for    phase    matching    between 
the    master    and     the    slave    beamformers.       Phase    matching    may   be 
obtained    for    an    arbitrary    array   geometry  by    augmenting    each 
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original  element  to  form  identical  subarrays  at  the  element 
locations.   Identical  preprocessors  may  then  be  used  to  form 
subarray  responses  with  nulls  in  the  selected  look 
direction.   The  preprocessor  outputs  must  be  cophased  for 
the  look  direction  and  applied  to  the  adaptive  beam  former 
while  the  original  element  outputs  are  cophased  and  applied 
to  the  slave  . 
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V.   SPATIAL  DITHER  ALGORITHMS 

A.   GENERAL  IDEA 

Spatial  dither  algorithms  [Ref.  2]  nave  been  newly 
conceived  for  the  purpose  of  applying  locally  controlled 
modulation  to  signals  arriving  at  angles  other  than  the 
Look-d i recti on ,  while  leaving  inputs  arriving  from  the  lo 
direction  unchanged.   The  technique  that  these  algoric 


re  employing  is  focused  on  a  reduction  of 


ca  amer    o  j  .,-  •; 


density    by    spreading     it    spectrally.       When    the    spati  2I     5it 
is    us  ed    with    a    conventional    ad  a  p  t  i  v  e   b earn  f  o nti er  ,     it    r  ed  uc  e  5 
the    signal    cancellation    effects. 

Conceptually,    a    simple    form    of    spatial    dither    algoricnm 
is    represented     in    Figure    5.1    and    stated    as    the    "3/4- inch 
plywood"    approach.       The    elements    of    an    antenna    array   may    oe 
imagined    to    oe    fixed    to    a    piece    of    plywood    which    provides 
rigid     insulating    support,    so    that    the    entire    array   may   be 
moved    mechanically.       In    either    one    or    two    dimensions,    the 
array    is    moved    in   directions    which    are    orthogonal    to    the 
look-direction.       Far-field    emanations    arriving    from    the 
look-direction    will    be    undistorted    by    the   mechanical    motion, 
while    emissions    from    off    axis    sources    will    be    distorted    by 
an    unusal     shi ft-of- time-base    form    of    modulation. 

The    outputs    of    the    antenna    elements    of    Figure    5.1    could 
be    applied    to    a    time    delay    and    sum    ( nonadapt  ive)     beamformer, 
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Figure    5.1.       Spatial    Dither    Algorithm     (3/4"    Plywood    Approach) 
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to  a  conventional  adaptive  beam  former,  or  to  a  Duval 1 
adaptive  beam  former .   Spatial  dither  could  be  beneficial  in 
each  case.   By  reducing  jammer  power  coherence,  so:ne  anti- 
jam  protection  is  provided  without  adaptive  beamforming,  an 
additional  anti-jam  protection  is  provided  with  adaptive 
beamforming.   Reduction  of  signal  cancellation  effect  in  a 
Frost  beam  former  can  oe  obtained  by  using  spatial  c.  i  ~ :".  _- ' 
preprocessing.   Breakup  of  jammer  signal  structure  is  3 
possiole  form  of  signal  preprocessing  applicaole  to  r.  1 '. 
types  of  adaptive  and  nonadaptive  beam  formers. 
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VI.   CONCLUSIONS 


^^  ^ 


.    IOa-; 


Adaptive  Noise  Cancelling  is  a  method  of  optimal 
filtering  that  can  be  applied  whenever  a  suitable  reference 
input  is  available.  The  principle  advantages  vf  the  met  ho.] 
that  an  A.N.C.  employs  are  its  adaptive  capability,  its  .  :\; 
output  noise,  and  its  low  signal  distortion.  T; 
capability  allows  the  processing  of  various  < 
properties  are  unknown.  Output  noise  and  signal  distort!  -i 
are  generally  lower  than  has  been  achieved  with  the  various 
conventional  optimal  filtering  configurations. 

The  FROST  beamformer  has  proven  itself  good  in  reject., 
various  interference  signals  and  enhancing  the  signal  of 
interest  for  several  cases.   interaction,  though,  between 
the  desired  signal  and  the  various  interference  signals  can 
lead  to  partial,  or  total,  cancellation  of  the  desired 
signal  within  the  adaptive  beam  former. 

The  DUVALL  composite  beamformer  gives  a  solution  to  the 
problem  of  signal  cancellation,  at  the  price  of  increased 
complexity  in  the  beamformer  structure  and  implementation. 
Under  that  scheme,  in  order  to  prevent  signal  cancellation, 
the  useful  signal  is  excluded  from  the  beamformer  in  which 
the  adaptive  process  now  takes  place  to  derive  a  set  of 
weights.   These  weights  are  now  copied  into  a  separate, 
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identical    "slaved"    processor    used    to    form    the    output    signal 
The    Duvall    approach    to    the    problem    is    appealing    but    since    i 
is    a    new   development,    possible    limitations    of    its 
performance    have    not    yet    been    assessed. 

A   WTDROW    solution    to    the    problem    involves    spatially,    or 
electronically,    moving    the    receiving    array    to    nod ul ate 
emanations    received    off    the    look-direction    without 
distorting    useful    signals    in    the    look-direction.       Tms 
approach,    called    "spatial    dither",     introduces    the    add:     .  )-  . 
possibility    of    modulating    or    smearing    "smart"    jammer 
signals,    thereby    limiting    their    effectiveness. 
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APPENDIX  A 
ANALYTICAL  SOLUTION  OF  THE  A.N.C. 

Consider  the  two-weight  noise  canceller  of  Figure  2.3, 
and  let  us  assume  that  the  samoled  reference  incuts  ire: 


xu  =  c  cos(vT 


: 


i  n  ( ui .  j  T  +  i> ) 


Figure  A. 1  is  a  flow  diagram  snowing  signal  propagati   :  i 
suc'n  an  adaptive  noise  canceller.   The  first  stec  to  '-   ie 
analysis  is  :o  obtain  the  isolated  impulse  response  fro;? 
error,  = ,  point  C,  to  the  filter  output,  point  G,  vith  th* 
feedback  loop  from  point  G  to  point  3  broken. 

Let  an  impulse  of  amplitude,  a,     to  be  applied  at  poirv 
at  discrete  time,  j  =  k;  that  is, 


e,  =  a«(j  -  k 


(A. 3 


The  6(j  -  k)  is  a  Kronecker  delta  function,  defined  as 


6(j  -  k)  = 


for  j  =  k 


otner wi se 


(A.  4 
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The    impulse    causes    a    response    at    point    D 


a   C    cos(m    kT    +    tO     for    j    =    k 


£JX1J 

ot'nerwi  se 


■,m:n    is    the    input    scaled     in    amplitude    by    the    instantaneous 
/alue    of    x,       at    j    =    k.       Tne    signal     flow    pa:.i    from    point    2    - 
point    E    is    that    of    3    digital     integrator     ,vith    transfer 
function,    2y/(z-l)  ,    and    impulse    response,    2uu(j-l)  ,    .vhere 
u  ( j )     is    tne    discrete    unit    seep    function 


u  v  ]  ) 


for    3    <    0 


otherwise  . 


Convolving    2uu(j-l)     with     e  T  *  tt    yields    a    response    at    point    E 
of 


W1J    =    2     M    aC    cos(11JokT    +    ^ 


(  A  .  7  ) 


where    ]>_-k    +    1.       When    the    scaled    and    delayed    step    function 
is    multiplied    by    x,,,     the    response    at    point    F    is    obtained 


y,        =     2     u     a    C       COS(a)    jT    +     |)cos(u    kT    +     t> ) 


(  A  .  3  ) 
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where    j    >    k    +    1.       The    corresponding    response    at    point    J, 
obtained     in    a    similar    manner,     is 


y2J    =    2    u    a    C       sin(jojT    +    ?)sin(^KT    +    p) 


( A .  3  ! 


where  j  >  k  +  1.   Combining  the  preceediny  two  equations 
yields  the  response  at  the  filter  output,  point  G: 


J  2: 


2     u  a  C   cos (  ,  T ( j-k) 


2  p  ^  C   u  (  j-k-1)  cos  ( to  T  ( j-k)  ) 


(A.  L 


Note    that    the    aoove    is    a    function    only    of     (j-k)     and    is    thus 
a    time    invariant    impulse    response    proportional    to    the     •     .    . 

impulse  . 

A    linear    transfer    function    for    the    noise    canceller    may 
now   be    derived    in    the    following    manner.       If    the    time,     <,     is 
set    equal    to    zero,    the    unit    impulse    response    of    the    linear 
time-invariant    signal-flow    path    from    point    C    to    point    G    is 


yj    =    2     n   C       u( j-1) cos( u0JT) 


( A .  1 1 


and  the  transfer  function  of  this  path  is 


G  (z)  =2  u  C" 


z  ( z  -  cos  a)  T) 
o 


z   -  2z  cos  uj  T  +  1 
o 


-1 
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2uC     (z    cos    Wqt    -    I 


A.  12) 


-    2z    cos    oj    T    +    L 


This    function    can    be    expressed    in    terms    o: 
frequency,     fl    =    2ir/T,     as 


a    r  ad  1  a; 


3d;ll  _>  . 


i  n  • 


g(z: 


2,  _i 

2  UC      [2     COS  (2  7TU      ft         )       -      1 


2  -1 

z       -     2  z    c  o  s  ( 2  it  -j    n       )     +    L 

o 


A .  I  3 


If  the  feedback  loop  from  point  G  to  point  3  is  new  cl 
the  transfer  function,  H ( z) ,  from  the  primary  input,  p 
A,  to  the  noise  canceller  output,  point  C,  can  oe  o'ota 
from  the  feedback  formula 


H(Z)     =    -— 


2z    cos ( 2  tw    2 
o 


_  i 


+    1 


~> 


2(1    -    iiC2)z    cos(2fj    n    l)     +    1    -    2uC" 


(  A .  1 4  ) 


This  equation  shows  that  the  noise  canceller  with  a  cosine 

reference  input  has  the  properties  of  a  notch  filter  at  the 

reference  frequency  along  the  signal  flow  path  from  primary 
input  to  output . 
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APPENDIX  B 
TAPPED  DELAY  LINE  FILTER 

The  tapped  delay  line  [Ref.  5]  is  the  preprocessor  most 
often  used  in  adaptive  filtering  applications.   Most 
mathematical  analyses  of  adaptive  filters  are  based  on  :•. 
tapped  delay  line  model.   A  capped  delay  line  can  be 


ie 


:    . 


implemented  either  in  analog  or  digital  form. 

of  an  adjustable  filter  based  on  an  analog  tapped  nelay  Line 

is  shown  in  Figure  3.1.   Such  a  delay  line  has  an  analog 
input  and  the  output  at  each  tap  is  a  delayed  version  of 

th 

this  input.   Thus  the  n  l  component  of  X(t)  is  the  output  at 

t'n        i     ■ 
tne  n ""  tap  ana  is  given  by 


xn(t)  =  x(t  -  (n  -  1)  A) 


(3.1) 


where  A  is  the  delay  between  taps  and  where 


X(t) 


x1(t) 

x2(t) 


xn(t) 


XN(t) 


(B.2) 
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,himm\    tlx«    daUy 


■  ft) 

^ 

***      X(t) 

T««Ot 

"l 

It) 

'J 

it) 

Wlqht    w 


H   £  HH   A 


s.,'*' 


Figure    B.l-       Tapped    Delay    Line    Filter 
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In  digital  implementation,  the  tapped  d.elay  line  is 
implemented  as  a  shift  register  either  in  hardware,  Figure 
B.2,  or  by  simulation  inside  the  computer,  Figure  3.3.   A 
new  digital  sample  of  the  filter  input  is  taken  at  every 
adaption  cycle  and  shifted  into  the  Leftmost  position  } f  t 
shift  register.   Old  samples  are  shifted  right  by  one 
position  and  the  oldest  sample  is  shifted  out  altogether. 
Each  position  of  the  shift  register  is  multiplied  by  the 
corresponding  weight.  '  Thus  the  digital  output  Df  t<vs  nil 
is  the  weighted  sum  of  the  last  N  samples  of  the  input. 
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APPENDIX    C 
COMPUTER    SOFISAPE    UTILIZED 


All    flcts    through    this    work      were    obtained    using 
and      C  S :'-!  P    language      programs    running      :n    the      I3.M    302 
f ram  e    ccnpu ter . 

The      following       CSM?      language      program      «as       use 
various    different   data)     to    produce    the    plots    2,5,     6, 
and    9.    As    inputs    to    the    progran    are    entered    the    coef. 
cf       the    transfer      function       H  (z)  =  Y  (zj  /X  (z)  in    the 

"-domain    ari   as      output    is    the    plot    of       the    magnitude 
frequency    cf    the    trarsfer   function    K(z). 


a  r  z  r  a  n 

r.  a  i  r. 

(  ~ '-  - 

1  2 

'    1  -I 

v  - :"  s  u  = 


//MANTKA1    JCE     (25o7,  1286)  ,  »  MANIKAS'  ,CLASS=E 
//*~CS:1A7    FE,DDN  AM  E=  t  IC7.  5Y  S VECTE, 0 ES T=10C A L 
//    Z"XEC    CSMFXV 
//X.SYSLI3    CD 
//  zc 

//  CD 

//  CD 

//  CD 

ESH  =  £YS3.,IMSI.5P,DISP=SHfi#VCL=SEE=MVS003rrJNIT=3350 

//X.CCMPEIN1    DD     EUMfll 

//X.  SSSFFINT    DD     DUMMX 

//X . EICIEAEH    DD    * 

SPICT    5CAIZ=.5     SEND 
//X.SYSIN    DE    * 

♦THIS    PRCGFAM    COMPUTES    THE     FREQUENCY    RL5?-y:iSZ    OF     UP 

*  TENTH    ORDER     FII1ER    WHOSE    TRANSFER     FUNCTION 

*  CAN    BE    EXPRESSED    AS; 


*** 


li- 


63 


*  H  (Z)  =  Y  (Z)/X  (2)  =  N(Z)/D  (Z)  , 

* 

*  -1  -2  -3  -u 

*  N(Z)=A0+A1*Z        +  A2*Z       +A3*Z       +  A4*Z       + 

* 

*  -1  -2  -3  -4 

*  D(Z)=EC  +  31*Z        +32*2       +33* Z       +  34*Z        + 


w  n  i  h  z 


+A10*Z 


10 
,ANE 


-  10 


+  3  1 0  *  Z 


HAGNI1UEE    AND    PHASE    ARE    GIVEN    FOR    THE    IKTE 
0     <    Yi    <    ?I       ONLY       SINCE     THAT    INFORMATION    I 
REPEATED    CN    THE    INTERVAL     PI    <    w    <2*?I. 
PtiASE     IS    CONS  TRAINED    TO     FLOS    OR    .UN  OS     U] 
THE    CFDINAIE     IS     It*    UNITS     C:     A/?I. 


5  V  Ai 


********************************************* 

INITIAL 

RENAME    IiaE  =  «OVRFI 

PARAM    EI  =  3.  141  59 

E  A  R  A  a    a  =  0 .  0  0 1 

PAR  AM    H  1=0, 25 

C1  =  -2.*CCS  (2.*PI*W1) 

C2=-2.*  (  1.-M)*COS  (2.  *EI*K1) 

C3=  1  .-2. *M 

A1=C1 

E1  =  C2 

E2=C3 

****  **********  *  *  **  *********  ****************** 

*  ENTER    THE     »A«    COEFFICIENTS 

PARAM    A0=  1.0,A2=  1.#A3=  0.,A4=  0 

PAR  A  a    A  6=  0.,A7=  0.,A8  =  0.,A9=  0 

**************  *************  ****************** 

*  ENTER     THE     '  E'     COEFFICIENTS. 

PARAM     E01.C  ,33=  0.  ,  B4=  0-,B5=  0 

PARAM    F6=  0.,37=  0.,E8=  0.,39  =  0 

********************************************* 


:**********  * 


***  *   * 

. ,  A  5  = 

.,A  10 
****  * 


.,3  10 
***  *  * 


******* 

* 

0  .* 

0.* 
******* 

* 

* 

0.* 
******* 
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EYNAKIC 

NOSCfil 

R=WOVEPI*PI 

EZGREE=1 80- *WO VEFI 
MAG 1=AO+A1*COS (W)+A2*COS [V* 
+A4*COS  (W*4.)  +A5*COS 
+A7*COS  (W*7.) +A8*COS 
+A9*COS (fl*9 .)  +A10*CO 
tv!AG2  =  A1  *SIN  (W)  +A  2*SI  N  (W*2.  ) 
+A  4*SI  N  (K * 4  . )  +  A5*SI N  (W 
+A7*SIN  (/T*7  .)  +A8*SIN  (W 
+  A9*SIN  (W*9  .)  +  A  10*SIN  ( 
MAG3=E0+E 1*C05  (W)+32*COS  (W* 
+  EU*CCS  (7J*4  .)  +E5*CCS 
+E7*CCS (W*7 . ) +E8*COS 
+E9*COS  (W*9.  )  +E19*CO 
HAG4=E1*SIN  (W)  +  3  2*SIN  (tf*2.) 
+3  4*SIN  (w*4.)  +B5*SIN  (W 
+  3  7 *S  I  N  ("7  *  7  . )  +  3  6  *  S I  N  { » 
+39*SIN  (»*9  -)  +3  10*SIS  ( 
MAG=S£ET  [MAG1**2+MAG2**2)/S 
IE  (MAG  1. EQ.0.0. AND.il AG2.EQ. 
IF  (MAG1.NE.0.0.OR.  MAG2.NE.  0 
ANG1=ATAN2  (M AG2 , M AG1 ) *  1 8 
IF(MAG3.EQ.0.0-AND.MAG4.EQ. 
IF  (HAG3.NE.0.0.  OH-  id  A  G  4 .  NE-  0 
ANG2=ATAN2 (MAG4/MAG3) *18 
EHASE=ANG1-ANG2 
IF (PHASE. GT. 180. )        PHASE=PH 
IF (PHASE. IE. -180 .)     PEASE=PH 
TIMEE    EINTIM=1.0 
PEIN1    W,EEG  BEE,  MAG,  PEASE 
CUTPUI    HCVRPI,KAG 
IABEI    DIGITAL       FREQUENCY       R 
IA3EI    CF    THE    NOTCH     FILTER     ( 


2. )  +A3*CG5  (W*3. )      .  •  , 

(U*5.)  +A6*COS  (W*6. ) 

( W  *  8  . ) 

S  (U*10.) 

+  A3*SIS  (W*3.)     ... 

*5- ) +A6*SIN  (S*6 . )     . . 

*8.)     ... 

N  *  1  0 .  ) 

2.)  +B3*COS  (W*3.)     .  . . 

(V'*5.)  +3 6* CO 5  (W*6. ) 

( rw  *  8  . ) 

S  ( K  *  1  0  .  ) 

+  B3*SIN(ff*3.J     ... 

*5-)+36*SIS(5r*6-J     .. 

*8.)     ... 

W  *  1  0 . ) 

CRT  (MAG3**2+MAG4**2] 

0.0)     AKG1=0. 

.  0)     ... 

O./PI 

0.0)     ANG2=0. 

.0)     ... 

O./FI 

ASE-360. 
ASE+360. 


ESPCNSE  -  MAGNITUDE 
WOVRPI=W/PI, 3  =  0.00  1) 
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EAGE     XYtlCT 

CUTEU1    KCVEEI, PHASE 

IA3EI    FREQUENCY       SESECNSE    -    PHASE      IS       DEGREES 

1A3EL    CE    THE    NOTCH    EIITEE     (  WOV5P  I=w  /PI,  'A-Q  .  00  1) 

PAGE    XYFICT 

END 

STOP 

ENDJCE 

/* 
/* 

/* 

/* 

The    following      FCBT&AH    program    was      used    to      produi 
tea  in  pattern    of    figure    3.5    . 

C 

C       THIS     EBCGBAM     DOSS    EC1AE    PLOTTING    Or     A    GIVEN    POLAE 

C       FUNCTION     (VTHE(I))     .VS.     THE    ANGLE     (3(1)) 

INTEGEE    I, J 

hIM    W  (73)  ,MAG  1  (73)  ,  H  AG2  (73)  , 

1  MAG 3  (73)  ,  KAG4  (73)  ,VTHE  (73)  , 

2  A0,  A1  ,A2,  £0,31,52,0  1,02,03  , 

3  MAXV/PIra#K1 
C 

PI  =  3.  14159 

M=.01 

W1=.125 

C1=-2.*C0S  (2-*EI*W1) 

C2=-2.*(1.-H)  *CCS  (2.  *EI*W1) 

C3=1. -2. *M 

A0  =  1. 

A1=C1 

A2=1  . 

E0  =  0- 

B1=C2 
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S2=C3 

EC    6    1    =    1 ,     73 

55(1)     =     (FICAT  (I)  -1.  0)  *5. 0*0.  0  17453 

6  CONTINUE 

£071=1,     73 

KAG1  (I)  =AC+A1*C0  S  (rw  (I)  )+A2*C0S  (W(I)*2.  ) 
MAG2  (I)  =A  1*SIN  (W  (I)  )  +A2*SIN  (W  (I)  *2  .) 
MAG3  (I)  =B0+3  1*COS  (:a  (I)  )  +  B2*C03  (W  (I)  *2.  ) 
MAG  4  (I)  =3  1*Si:,T  (W  {I)  )  +32*  SIN  (H  (I)  *2.) 

7  CONTINUE 

EC     10    I    =    1,     7  3 

VTH  3  ( I)  =  S  Qvl  { HAG  1  (I )  *  *2+  MA  32  (I)  *  *  2 )  / 
5  S  GET  (MAG  3  (I )  *  *  2+  MA  ^4  ( I )  *  *  2 ) 

V:(l)     =    W  (1)     *    57  .29  6 
10  CONTINUE 

C    ***    NCEMAIIZE     AND     BCHM    THE    CDTPUI    ARRAY    VTHZ(I)      *** 
C 

MAX7=VTKE(1) 

EC     12    I     =    2,     73 

IE     (  (VTHE  (I)  )  -LE.  (MAX?)  )     GO    TC    12 
MAXV  =  VTKE  (1) 

12  CONTINUE 

DC    13    I    =    1,    73 

VTHE  (I)     =    TITHE  (I)  /WAXV 

13  CONTINUE 
C    SET    UP    TEK618 

CAII    TEK618 
C  CAIL   COMPRS 

CAII    PAGE(8. 5,11.) 

CAII    BIO  WUP  (.6) 

CAII    N03SDR 

CAII   A  RE A2D  (7.5,7.5) 

CAII    HEAEIN     ('BROADSIDE    ARRAY    P ATTEEN ' ,- 2 3, 2 . ,  1 
C    SET    UP    POIAR    GRID 

CAII    POIAR  (0.01745,. 26 6,3. 75, 3. 75) 
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CAI1  MAfiKSB ( 16) 

CALL    DOT 

CALL   CUBVE  (W  ,VTHE,73,  1 

CALL    RESIT      (  'EOT') 

CALL    3EID(2  ,  4) 

CALL    ENDPL (0 ) 

STOP 

END 
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